A simple radionuclide-driven single-ion source 
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We describe a source capable of producing single barium ions through nuclear recoils in radioactive decay. 
The source is fabricated by electroplating 148 Gd onto a silicon a-particle detector and vapor depositing a 
layer of BaF2 over it. 144 Sm recoils from the alpha decay of 148 Gd are used to dislodge Ba + ions from the 
BaF2 layer and emit them in the surrounding environment. The simultaneous detection of an a particle in 
the substrate detector allows for tagging of the nuclear decay and of the Ba + emission. The source is simple, 
durable, and can be manipulated and used in different environments. We discuss the fabrication process, 
which can be easily adapted to emit most other chemical species, and the performance of the source. 

PACS numbers: 29.25.Ni, 34.50.-s, 81.15.Pq 



I. INTRODUCTION 

It is often necessary in experimental settings to pro- 
duce specific ionic species in a controlled manner. In- 
deed a large number of techniques to produce ions have 
been documented, from the early Penning sources to laser 
and electron-beam produced plasmas to thermal ioniza- 
tion^. These sources, however, can only operate under 
vacuum, and in many cases require expensive and com- 
plex infrastructure to allow for ion separation and deliv- 
ery. 

The work discussed here arises from the need to de- 
velop a simple source of monoatomic, singly-ionized bar- 
ium that can be used in a cryogenic liquid, a gaseous 
environment, and in vacuum. Such a source is necessary 
for the R&D efforts of the Enriched Xenon Observatory 
(EXO). The EXO Collaboration is planning a series of 
experiments designed to determine the mass of the neu- 
trino by searching for the neutrinoless double-beta decay 3 
of 136 Xe. To suppress possible radioactive backgrounds, 
EXO plans to add to the standard low background tech- 
niques that are customary in these experiments the abil- 



a ) Now at Indiana University, Bloomington IN, USA 

b ) Now at the University of Seoul, Seoul, Korea 



ity to tag the barium daughter of the double-beta de- 
cay using resonant fluorescence 4 ^. Since possible detec- 
tor technologies would include xenon Time Projection 
Chambers (TPCs) in high pressure gas^or liquicP phase, 
the R&D work on Ba-tagging requires a method of releas- 
ing Ba + in such media. The injection through a thin win- 
dow of Ba ions produced in an accelerator would clearly 
be possible but it would require expensive equipment and 
would result in very highly ionized states, owing to the 
electron-stripping in the entrance window and in the high 
density medium. 

The source presented here makes use of the Sm daugh- 
ter of the decay 

64 8 Gd-> JfSm + a (1) 

to dislodge a Ba + ion from a BaF 2 layer coated over 
the Gd source. The a recoiling against the Sm is then 
detected in a semiconductor detector that is used as a 
substrate for the assembly, providing a tag for the ion 
emission. 148 Gd is a convenient a emitter because its 
half-life (74.6 y) provides a good compromise between 
high specific activity and durability. In addition, the 
decay of 148 Gd is directly in the ground state of 144 Sm, 
so that no other radioactivity is produced. Finally, Gd 
can be reliably electroplated in thin, uniform layers, as 
will be discussed. 
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The Gd source is deposited onto a PIPS® a detec- 
tor^ and a special technique is used to confine the ac- 
tivity to a small central region. A thin layer of BaF2 
is then vapor-deposited over the entire detector. The as- 
sembly is mounted in a custom metal-ceramic holder that 
is suitable for ultra-clean operations and moderate vac- 
uum bakeouts. The ion source obtained is inexpensive, 
compact and suitable for a variety of environments. The 
selection of coatings other than BaF2 can provide similar 
sources for most elements. 



II. SOURCE FABRICATION 

148 Gd is obtained in the form of GdCl3 in HC1 aque- 
ous solution. Natural Gd is present as carrier in the solu- 
tion in the approximate ratio 9 nat Gd: 148 Gd=l:6.6x 10" 4 . 
Gd electroplating works best in isopropanol (IPA) from 
gadolinium nitrate^. So the chloride is transformed into 
nitrate as: 

GdCl 3(aq) + 3HN0 3 ( aq ) -> Gd(N0 3 ) 3 (a q ) + 3HCl (aq) 

(2) 

The reaction is achieved by drying a 500/iL (0.1/iCi) 
148 GdCl3 batch on a hot plate and adding 500/iL of 1M 
HNO3 11 . The water is then evaporated again and more 
HNO3 is added, repeating the process three times. Af- 
ter the last drying cycle, the 148 Gd(N03)3 is dissolved in 
500/iL of IPA 12 . 25/iL of this solution are calculated to 
have an activity of 185 Bq. A measurement of the activ- 
ity, performed by dissolving one 25/iL batch in a liquid 
scintillation counter 13 , confirms a carry-over efficiency of 
~100%. Since the Gd(N03)3 concentration is extremely 
low, a drop of 0.05M HN0 3 is added in 10L of the IPA in 
order to increase the conductivity of the plating solution 
to an acceptable level. 

The use of a conventional plating cell is excluded be- 
cause of the requirement to deposit the 148 Gd only in a 
specific and small region of the detector. Therefore, a 
plating circuit is established by mounting a PIPS® de- 
tector on a horizontal support and lowering an anode to 
a position about 0.5 mm above the center of the detector 
using a micrometer, as shown in Figure [I] Using a cali- 
brated pipette a small drop of ~ 5/iL of solution is then 
deposited between the anode and the detector, which 
serves as the cathode in the circuit. The drop remains 
in place because of surface tension. A 200 V potential is 
then applied to the circuit via a 10 current-limiting 
resistor, establishing a current that is measured to be 
18 /iA. Depending on the activity required more 5/iL 
batches are added, as the solvent evaporates, for a total 
of 25 to 50/iL. Finally, pure IPA aliquots are added, again 
to counter evaporation. A total plating time of 15 min 
was found to be adequate by trial and error by measuring 
the activity on the detector (since the charge transport 
in the cell is completely dominated by the HNO3 added 
to the IPA and has no relation with the amount of Gd 
deposited). Because of the non-radioactive Gd carrier, a 




FIG. 1. Schematic view of the radioactive Gd plating setup. 
The micrometers allow for an accurate and stable positioning 
of the anode ~ 0.5 mm above the center of the detector being 
plated. A drawing of the Canberra PIPS® detector used is 
provided in the inset. 



source of ~200 Bq activity, obtained from 25/iL plating 
solution, roughly covers a surface of 8 x 10 -2 cm 2 with 
~10 monolayers of Gd. 

Initial tests using Pt anodes were found to produce 
large contaminations in the Gd coatings. This was evi- 
dent by optical inspection of the plated detector. Such 
contaminations are unacceptable because they reduce the 
energy of the Sm recoils, as confirmed by the analysis of 
the a energy spectrum from the source measured in an 
external surface barrier detector. Anodes obtained using 
4 mm - wide slices of a B-doped silicon wafer (resistivity 
« 50 fhn) were found to be sufficiently clean and were 
used for the work presented here. 

The a detector used as a substrate has surface con- 
tacts formed by shallow ion implantation and it is espe- 
cially appropriate because of its ruggedness. Its hexag- 
onal shape (also shown in Figure [T]) is of no particular 
significance here. The functionality of the detector af- 
ter plating and the quality of the plated source are ver- 
ified by a-counting both with the source-substrate de- 
tector ("self-counting") and, in some cases, with an ex- 
ternal surface-barrier detectoi^ in a vacuum chamber. 
In addition, an external 241 Am source is used to verify 
the detector performance and resolution after plating. A 
typical 241 Am spectrum recorded by the PIPS® detec- 
tor with a 50 V bias and 2 /is shaping time, before and 
after Gd plating is shown in Figure [2] The three peaks 
at 5485.6 keV, 5442.8 keV and 5388.2 keV are clearly 
visible before plating, with a 22 keV FWHM resolution 
for the most prominent peak. This is somewhat worse 
than the 14 keV FWHM stated by the manufacturer. 
The spectrum after Gd plating is substantially broader 
(~ 130 keV FWHM), presumably due to some deterio- 
ration of the detector and some energy loss in the Gd 
and other impurities left on the surface by the plating 
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FIG. 2. a spectrum of an external 241 Am source recorded by 
the PIPS® detector before (full circles) and after (empty cir- 
cles) the Gd plating. The three main a peaks at 5485.6 keV, 
5442.8 keV and 5388.2 keV are clearly distinguishable in the 
first data set. In this case the resolution of the 5485.6 keV 
peak is of 22 keV FWHM. After plating, the energy resolution 
degrades substantially and the new resolution of ~ 130 keV 
FWHM makes it impossible to resolve the three 241 Am lines. 



process. This degradation, however, appears to be stable 
in time and is not important for the simple a detection 
required in tagging the ion emission. The spectrum from 
the single 3182.7 keV a emitted by the 148 Gd plated on 
a detector shows (Figure [3^a)) a resolution of ~ 60 keV 
FWHM. 

The activity of the plated 148 Gd is measured by read- 
ing out the signal from the a detector, as will be de- 
scribed in the next section. For typical sources this ac- 
tivity is ^200 Bq and matches, within experimental un- 
certainties, the activity measured in the plating solution, 
implying that the plating efficiency is rather high. 

A coating of the appropriate chemical species and 
thickness will generally produce free ions and neutral 
atoms dislodged by the collision of the 89 keV recoiling 
144 Sm. This is expected to be true for most elements, 
although with varying efficiencies, mainly due to the ra- 
tio of masses between the Sm and the "target" species. 
Different (in particular, heavier) a emitters (e.g. 208 Po 
or 209 Po) could replace the 148 Gd, although this was not 
attempted in our work. 

BaF2 was chosen as the chemical species containing 
Ba because of its relatively good stability in aii^ and 
the fact that the other atomic species, fluorine, has a 
mass that is very different from that of Ba. The optimal 
layer thickness is estimated using the SRIM 16 simulation 
package with full damage cascades simulation. Figure [4] 
shows the predicted Ba yields as a function of thickness. 
BaF 2 was simulated by simply providing SRIM with the 
appropriate number densities of Ba and F. Yields greater 
than one, like those reported in Figure [4j are not unusual 
in sputtering processes. 

While the Ba yield can be larger than 1, as in the 
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FIG. 3. 148 Gd "self-counted" a spectra. Panel (a) shows 
the spectrum of the 148 Gd layer only. The FWHM energy 
resolution is ^60 keV. Panel (b) shows the spectrum after 
40 nm of BaF2 are deposited over the source. The FWHM 
energy resolution is in this case ^170 keV. 



case of sputtering, it should be kept in mind that SRIM 
does not distinguish between neutral atoms and ions in 
different charge states, so that the yield of Ba + is not 
known from the simulation. The energy spectrum for 
Ba is predicted by SRIM to be peaked at zero with a 
long tail extending above 100 eV. 15 nm was chosen as 
the optimal layer thickness for the BaF2, although some 
early sources were fabricated with thicknesses from 8 nm 
to 40 nm. 

As already mentioned, electroplating was chosen as the 
technique to deposit the a emitter because of the layer 
quality and the possibility of obtaining deposits that are 
easily limited to a specific area. The technology for the 
overcoating by the target species was then chosen so as 
not to chemically interfere with the radionuclide and, at 
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FIG. 4. Yield for Ba atoms as a function of layer thicknesses, 
as calculated by a SRIM simulation. Since SRIM does not 
distinguish between different charge states the yield from the 
graph includes neutral atoms, singly charged ions as well as 
ions with a higher ionization state. 

the same time, provide a uniform layer over the entire 
front surface of the detector. Vacuum evaporation is suit- 
able for this since the energy of the ions is low (compared 
to sputtering) and no wet chemistry is involved. In our 
system, tantalum boats are ohmically heated to an ap- 
propriate temperature and the thickness of the deposited 
layer is monitored by a microbalance located next to the 
source being plated. Most materials can be deposited 
in this way or using electron-beam heating for refrac- 
tory materials. As shown in Figure [3^b), the a energy 
resolution from the 148 Gd source is further deteriorated 
after the BaF2 coating. While the origin of the deterio- 
ration is not understood, the detector is still functional 
and completely adequate to tag the decay and the sub- 
sequent emission. In order to ensure the best possible 
purity for the BaF 2 coating, a small scintillation-grade 
BaF 2 crystal was used in the evaporation process. 

Sources prepared as described above are then mounted 
in the special holders shown in Figure [5] All metallic 
parts are made of gold-plated stainless steel, and the insu- 
lating ring separating the two diode contacts is made out 
of Macor®. The elastomer used in the standard encap- 
sulation to maintain a constant pressure on the contacts 
is replaced by a gold-plated wavy washer. The holders 
allow for cooldown to 77 K and are compatible with the 
extreme cleanliness required for ultra-high vacuum oper- 
ations as well as for operation in liquefied noble gases (in 
particular liquid xenon). 

III. SOURCE CHARACTERIZATION 

The ion emission from the source is analyzed with a 
simple time of flight (TOF) spectrometer comprising two 
Einzel lenses and designed to maximize the transport ef- 
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FIG. 5. Top: Exploded view of an ion source in its metal- 
ceramic UHV mount. The "wavy washer" provides a uniform 
force on the detector to ensure proper contact. All metallic 
parts are made out of gold-plated stainless steel. The ring 
to the right of the drawing is made out of Macor® . Bottom: 
Photograph of a completed source. 



ficiency of ions from the rather large emittance produced 
by the source to a channel electron multiplier^ (CEM). 
The spectrometer is designed with a central drift region 
and electrodes at the two ends that allow for the indepen- 
dent adjustment of the electric field at the surface of the 
source (Source, 350 V/cm < E source < 750 V/cm) and 
the ion impact energy on the CEM (K 2400 eV < K < 
3800 eV). The overall length of ~ 17 cm produces a TOF 
for Ba + of ~5 /xs, depending on the exact configuration 
of potentials. The "start" signal for the TOF measure- 
ment is provided by the PIPS® detector on which the 
source is built and the "stop" is provided by the CEM 
pulse. In this case the 45 V bias of the PIPS® detector 
is provided by a battery, so that the entire source can 
be floated to high voltage as required by the different 
configurations of the spectrometer. Both signals from 
the source and the CEM are digitized with flash ADCs 
(FADCs) at 250 MS/s. The signal from the source is 
used to trigger an acquisition. While a fast preamplifier 
is used for the CEM pulse, the signal from the source 
is fed into an Ortec model 142 preamplifier followed by 
a model 474 timing filter amplifier (TFA) with 100 ns 
shaping. Timing thresholds are applied offline. 

A typical TOF spectrum obtained for the intermedi- 
ate values of Source = 550 V/cm and K = 2800 eV is 
shown in Figure [6| The prominent peak around 5 /is has 
a position and width consistent with those predicted by 
a SIMION 19 simulation for Ba + ions. Such simulation 
includes the correct natural isotopic mix for the atomic 
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FIG. 6. Time-of-flight spectrum for a typical Ba source. The 
source has a measured 148 Gd activity of 173 Bq and a 15 nm 
thick BaF2 coating. The Ba + yield, measured as the integral 
of the main peak, is of 2.7 Hz. 



mass of Ba and an initial energy distribution derived by 
the SRIM output, as explained below. A fit of this peak 
to five Gaussian functions representing the main Ba iso- 
tope masses provides an estimate for the spectrometer 
FWHM resolution of 5.8 amu at ^137 amu. The assign- 
ment of the peak to Ba + is corroborated by the observa- 
tion that no such peak is present in a source where the 
BaF2 coating is replaced by Al. In addition, and maybe 
most importantly, Ba deposited from a source of this type 
onto a clean silicon substrate was then desorbed and res- 
onantly ionized with two lasers. From this very specific 
result it is clear that positive Ba ions are indeed emitted 
by the source. 

According to SRIM the energy spectrum of the Ba 
emitted by the source peaks at zero and decreases mono- 
tonically, with 17% of the spectrum having energies 
greater than 300 eV. Since, as already mentioned, the 
SRIM simulation is unable to distinguish between neu- 
tral and charged states, a different method is required 
to estimate the expected ion fraction, a = n + /n to tai- 
Ideally this ratio would be known as a function of Ba en- 
ergy. Unfortunately, a general understanding of the ion 
fraction in sputtering is still an open question, even for 
pure metals, and data is completely lacking for ionic com- 
pounds, such as BaFj-^4 Some models and experimental 
data describe a by a decaying exponential in the inverse 
velocity plus a constant term, with the transition tak- 
ing place around a critical velocity of ~2 cm//is 20 . The 
net result of this function combined with the sputtering 
energy distribution is that the fraction of ions is low at 
low energy and approaches unity at higher energy. The 
transition takes place around 2 cm/ /is, or ~300 eV in the 
case of Ba. 

The actual energy distribution of Ba + ions affects 
the transport efficiency of the spectrometer significantly. 
While such efficiency is calculated to be 85% for the en- 



tire energy distribution produced by SRIM, it is only 8% 
for the part of such energy distribution exceeding 300 eV. 
Also the TOF peak position and width are to some extent 
predicted to vary for different regions of the emittance of 
the source. A simple consistency check for the Ba + yield 
can be made by assuming in the simulation that only ions 
(neutrals) are produced above (below) a certain energy 
cutoff. The yield, TOF peak position and width are then 
compared with data for a variety of values of Source and 
K. A reasonably good description of the data occurs for 
a 300 eV cutoff. The simulated rate reproduces the data 
within a factor of two when E source and K are varied from 
550 V/cm and 2800 eV to the extremes mentioned above. 
With the same cutoff value the simulation also provides 
a good prediction for the Ba + TOF peak position (bet- 
ter than 10%) and width (better than 25%), again when 
scanning Source an d K from their central values to the 
extremes. 

For the spectrum in Figure [6] the Ba + yield is measured 
to be 2.7 Hz. Using the a activity of 173 Bq, SRIM 
predicts a rate of Ba emission of 430 Hz for the 15 nm 
thick BaF2 coating (see Figure |4|. Assuming that the 
ions comprise only those Ba atoms ejected with energy 
above 300 eV, or, according to SRIM 17% of the total, the 
Ba + emission rate is calculated to be 73 Hz. Using the 
spectrometer efficiency of 8% appropriate for this energy 
range, a rate of 5.8 Hz is predicted to arrive at the CEM. 
The ratio between this figure and the measured rate of 
2.7 Hz is consistent with typical CEM efficiencies for ions 
with the velocities expected here. 

TABLE I. Possible assignment of the TOF peaks for the 
source in Fig. [6] The associations of H + and Ba + (in ital- 
ics) are assumed and used to compute the mass values for the 
other peaks. 



Time 


Measured mass 


Possible assignment 


(fJS) 


(amu) 


(known mass in [ ]) 


0.31 


1.01 


H + [1.01] 


1.44 


12.9 


C + [12.0] 


1.99 


23.5 


Na + [23.0] 


2.16 


27.6 


Si + [28.1] 


2.23 


29.3 




2.75 


43.5 


C 2 H 5 + [45.1] 


3.19 


57.9 


C 3 H 8 + [60.1] 


3.55 


70.9 


Ba 2+ [68.7] 


4-99 


137.3 


Ba + [137.3] 


5.33 


156.2 


BaF + [156.3] 



The TOF spectrum for light masses is similar to the 
one observed for the source coated with Al. Its interpre- 
tation nevertheless poses some challenges because of the 
limited resolution of the spectrometer, the large emit- 
tance of the source and the difficulty in predicting the 
charge to neutral ratio as a function of the species and 
initial energy. A possible assignment, obtained using as 
reference the Ba + peak and the peak at 0.31 /is, assum- 
ing to be due to H + is shown in Table [i] The presence of a 
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prominent Na + peak is not surprising, given the common 
presence of sodium contamination on surfaces. Emission 
of some Si from the substrate is not unexpected, while 
C2H50 + and CsH80 + are commonly observed in mass 
spectra of the IPA that is used in the Gd plating pro- 
cess. No compelling assignment is available for the peak 
at 2.23 /is. Finally, the satellite peak at a mass larger 
than that of Ba is not present in the Al coated source 
and it may be attributed to BaF + . 



IV. CONCLUSIONS 

We have built a simple, inexpensive and compact 
source of single Ba + ions using the recoils of an a emit- 
ting radionuclide to dislodge Ba atoms from a BaF2 coat- 
ing. The source uses a surface barrier detector as a sub- 
strate, so that a tag is provided for the ion emission. Its 
simplicity and the possibility of delivering tagged ions 
in various environments make it attractive in applica- 
tions where a low rate of individual ions is desirable. In 
addition the emission of neutral Ba, not useful for our 
purposes, may find other applications. The replacement 
of the BaF2 with other materials provides the means of 
emitting ions from a large number of different species. 
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